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ABSTRACT 

We study semi-analytically and in a consistent manner, the generation of a mean 
velocity field U by helical MHD turbulence, and the effect that this field can have on 
a Mean Field Dynamo. Assuming a prescribed, maximally helical small scale velocity 
field, we show that large scale flows can be generated in MHD turbulent flows, via 
small scale Lorentz force. These flows back-react on the mean electromotive force 
of a Mean Field Dynamo through new terms, leaving the original a and (3 terms 
explicitly unmodified. Cross-helicity plays the key role in interconnecting all the effects. 
In the minimal r closure that we chose to work with, the effects are stronger for large 
relaxation times. 
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1 INTRODUCTION 

Magnetohydrodynamical (MHD) turbulence seems to be a 
major physical process to generate and maintain the mag- 
netic fields observed in most of the structures of the Uni- 
verse (Branden burg fc Subramania n 2005a; Zel'dovich et al.l 
1 19831 ). When addressing the problem of the generation of 
large scale magnetic fields by small scale turbulent flows, 
a model kn own as Mean Field Dynamo (MFD) is usually 
considered (|Moffatj[l97ah . Despite its simplicity and lack 
of broad applicability, it proved to be a very useful tool in 
studying qualitatively conceptual issues of large scale mag- 
netic field generation. The mechanism is based on decom- 
posing the fields into large scale, or mean fields, U, B, 
A and small scale, turbulent ones u, b, a. These small 
scale fields have very small coherence length, but their in- 
tensities can be higher than the one of the mean fields. 
In this mechanism the evolution equation for B is written 
dB/dt = V x (UxB + £- t)j), where J = V x B, rj is 
the Ohmic resistivity and £ = u x b the turbulent electro- 
motive force (t.e.m.f.) 0. The term U x B is usually disre- 
garded in the studies of MFD as the focus of most of them 
is to understand the generation of large scale quantities due 
to small scale effects. If homogeneous and isotropic turbu- 
lence is considered, the t.e.m.f can be written as £ — aB — 
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1 Overlines denote local spatial averages: they represent vector 
quantities whose intensities may vary in space, but whose direc- 
tion and sense are uniform or vary smoothly. ()„„i denote volume 
averages, i.e., quantities that can depend only on time. 



/3J, with a 



- (1/3) Tcorr [u ■ (V x u) - b • (V x b)J and 



1/3) TeorrM 2 , Tcorr being a correlation time (iMoffatt 



| l972j : iRudiged 1 1974 |Pouquet_at "ail 1 19761 ; IZel'dovich eta! 



Il983f ). The dependencies on B and b are due to t he back- 
react i on of those induced fields on the dynamo ((Moffatt 
ll97StlBrandenburg fc Subramanianll2005al ). In the kinemat- 
ically driven dynamo considered here, the features of the 
generated fields crucially depend on the helicity of the flows: 
helical flows are at the base of the mechanisms to generate 
large scale fields, while non-helical flows would only produce 
small scale fields. This separation, however, is somewhat ar- 
tificial, a s small scale fields are also produced by helical tur- 
bulence (jBrandenburg fc Subramanianll2005al ). 

In this paper we want to address an issue not (or very 
seldom) considered in the literature up to now, namely, 
the induction of large scale flows U, also named shear 
flows, by the small scale turbulent fields, and how these in- 
duced flows back-react on the turbulent electromotive force 
£ of a MFD. On one side, we mean that the expression 
dB/dt — V x (£ — r)3} would be valid only during the time 
interval in which U x B <C £\ and on the other, even if 
this conditions is satisfied, £ could be affected by the gen- 
eration of U and consequently its functional form should 
be modified to incorporate this effect. The generation of 
magnetic fields due to the action of these large scale ve- 
locity flows instead of by £ was recently addressed ana- 
lytically by several authors JRogachevskii fc KleeorinI 120031 . 
12004 Radlcr fc Stepanovl 2005). and was also studied nu- 
merically by iBrandenburgl (|200l|) and semi-analytically by 
IBlackman fc Brandenburg (|2002l ). However, none of those 
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works addressed specifically the issue we want to analyze 
here. 

We work in the framework of the two scale approxi- 
mation, that consists in assuming that mean fields peak 
at a scale k^ 1 while turbulent ones do so at fc^ 1 <C k^ 1 , 
and also consider homogenous and isotropic turbulence. Al- 
though this kind of turbulence is of dubious validity when 
dealing with large scale fields, it serves well for initial, qual- 
itative studies of the sought effects. Another assumption we 
shall make is that B is force-free, i.e., of maximal current 
helicity. Although fields with this feature can be observed in 
certain astrophysical environments, they are not a general- 
ity, and also they are not seen in some numerical simulations. 
The main reason to use them here is to simplify the (heavy) 
mathematics, while maintaining a physically meaning sce- 
nario. 

In order to find £ when Lorentz force acts on the 
plasma, we must solve a differential equation that contains 
terms with one point triple correlations, i.e., averages of 
products of three stochastic fields evaluated at the same 
point. This means that instead of dealing with only one 
equation to solve for £, we have to solve a hierarchy of them. 
In order to break this hierarchy and thus simplify the math- 
ematical treatment of the problem we must choose a clo- 
sure prescription, which consists in writing the high order 
correlations as functions of the lower order ones, but main- 
taining the physical features of the problem under study. 
In MHD the intensity of the non-linearities is measured 
by the magnetic Reynolds number, which is defined from 
the induction equation for the magnetic field, dJi/dt = 
V(U x B) + r?V 2 B, as R m ~ (UB/l) / {r)B/l 2 ) = Ul/f], 
where now U is a characteristic velocity of the plasma, 
I is a characteristic length and r\ is the ohmic resistivity. 
Thus for R m <C 1 the non-linear terms can be neglected in 
front of the resistive ones, and therefore the equations be- 
come linear, while for R m 3> 1 the resistive terms should 
be dropped off in front of the non-linear ones. The inter- 
mediate regime is more difficult to analyze. In this paper 
we shall consider R m 3> 1, i.e., the non-linearities must be 
maintained, and consequently a closure scheme must be se- 
lected to deal with them. We ch oose to work with the so 
called minimal r approximat i on ([Blackman fe F ield 2002; 
iBrandenburg fc Subramanianl l2005al ). whereby the triple 
moments in the equation for £ will be considered as pro- 
portional to quadratic moments, and written in the form 
(,£, with the proportionality factor f ~ T rel' where r re ; 
is a relaxation time that can in principle be scale, and/or 
R m , depend ent. The validity of this closure was checked 
numerically l|Brandenburg fc Subrama nian 2005b) for low 
Reynolds numbers, and was also verified for the case of pas- 
sive scalar diffusion. We shall assume that it is also valid 
for all the triple correlations that may appear throughout 
this work. We consider boundary conditions such that all 
total divergencies vanish. These conditions may be a bit un- 
realistic for astrophysical systems, but they have two advan- 
tages: magnetic helici ty becomes a gauge-invariant quantity 
jBerger fc Fieldlll984h . and the obtained results can be com- 
pared with numerical simulation, as to perform them it is 
customary to use those conditions. We consider fully helical, 
prescribed, u fields. 

Our starting points are the evolution equations for £ 
and for the magnetic helicities H^s, as the evolution of 



these quantities is tightly interlinked (|Blackman fc Field! 



2002; iBrandenburg fc Subramanianl |2005a[ ) in the absence 
of shear flows. When including U new terms appear in 
the equation for £ , but the ones that drive the e volution 



of £ in absence of U, i.e., the a and terms (IM offatt 



1972; Rudigcr 1974; iPouquet at al. Hl976l ; lBlackman fc Field! 



2002J), are not explicitly modified. For those new terms, fur- 



ther equations must be derived, that in turn show the sub- 
tleties of the interplay among U, u, b and B. Due to the 
chosen boundary conditions, the evolution equations for Hjf 
and H^ 1 do not explicitly depend on U, they will do so 
implicitly through £. We consider fully helical U fields, so 
we study their growth through its associated kinetic helic- 
ity, H u = ((V x U) -U) and show that, for fully he- 
lical, prescribed u, large scale flows will be always gener- 
ated, as long as small scale Lorentz force is not null, i.e., if 
(Vxb)xb^O. We shall consider two values for the mag- 
netic Reynolds number, R m = 200 and 2000, and for each 
case analyze the effect of short and large r re ;. In general we 
find that for short r re i (large Q, i.e., strong non-linearities, 
the effect of large scale flow is negligible, thus producing 
results that practically do not differ from the ones in the 
absence of large scale flows. For large r re i (small f) the gen- 
eral effect is an enhancement of the electromotive force and 
the inverse cascade of magnetic helicity, this enhancement 
being stronger for R m = 2000 than for R m = 200. 



2 MAIN EQUATIONS 

Ohm's law for an electrically conducting fluid reads E = 
— U x Ji + r/J, with 77 the electric resistivity and J the electric 
current. The equation for B is the induction equation: 



<9B 



dt 

and from B 
E= -V$- 



V x E, 



= V x A we have 
dA 
~dt' 



(1) 



(2) 



which is an evolution equation for A. The equation for the 
velocity field U is the Navier-Stokes equation that, when 
considering only Lorentz force, reads 

= - (U • V)U-^E + (V x B)xB-/A7x(V x U).(3) 

with v the kinetic v iscosity. To work within mean field 
theory (Moffatt 1978J) we decompose the different fields as 
B = B+b, A = A+a, U = U+u, E = E+e and * = ¥+0, 
where any mean value of stochastic quantities vanishes. 
The derivation of the evolution equations for the mean and 
stochastic fields i s a standard procedure, already described 
in th e literature (|Zel'dovich et al.l 1 19831 ; iBlackman fc Field] 
2002). Consequently we only write here the results. Assum- 
ing incompressibility of the large and small scale flows, con- 
sidering that B is force free and working with Coulomb 
gauge for the vector potential, i.e., V • A = = V ■ a, 
we obtain the following equations for the mean fields 



<9B 



^ = V x rUxB + £-r, VxBl 
Ot J 

^=PPtJxB + £l-T,VxB, 
ot 1 J 



(4) 
(5) 
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where £ = u x b is the t.e.m.f., and 
~dt 



P [- (U ■ V) U - (u ■ V) u + (b • V) bj 



+ zATU. 



(6) 



(p) . . = Sijd 2 — didj is the projector that selects the sub- 
space of solutions of eq. ([2]) that satisfy the Coulomb gauge 
condition and the subspace of solutions of (0 ) that satisfy 
the incompressibility condition. Observe that eq. ((6} shows 
that a large scale velocity field can be induced from an ini- 
tially zero value, as long as — (u ■ V) u+ (b ■ V) b ^ 0. The 
equations for the small scale fields read 

^ = V x |U x b + u x B + u x b -I] +rjX7 2 h, (7) 
dt L J 

9a 



— = P[Uxb + uxB + uxb-J] + r/V 2 a, 



dt 
and 
du 



(8) 



^— = P I- (U- V) u (u- V)U- (u- V)u+ (u- V)u 

dt 



+ (B ■ V) b + (b ■ V) B + (b ■ V) b - (b ■ V) bj 



+ uV u. 



(9) 



2.1 



Evolution Equation for Derived Quantities: 
Magnetic Helicity, Large Scale Kinetic 
Helicites, and the Stochastic Electromotive 
Force 

As stated in the Introduction, we want to study if an initially 
zero, or very weak U, can grow due the action of a MFD, 
and back-react on it and on the magnetic helicity. This last 
quantity is defined as the average over t he entire volume 
of the dot product A ■ B (Biskamp 1997). In this way, we 
write the magnetic helicity associated to the large and small 
scale fields respectively as 1 = (A • B^ ( and Hg 1 = 
(a ■ h) vol , and by definition they can only depend on time. 
The evolution equations fo r Hf 1 and , for the chosen 
boundary conditions, read dBlackman fc Fie ld 2002) 



dt 
and 



dt 



2 <£- B L i -M(VxB).B) uo; , 



= - 2 ( £ - B )^- 2 ^( Vxb )- b )^ 



(10) 



(11) 



Observe that these equations have the same form as the ones 
obtained in the absence of large scale flows. This fact is due 
to the selected boundary conditions: magnetic helicity can 
be injected into the system through the boundaries by large 
scale flows. Thus, in the case under consideration here, these 
flows cannot explicitly transport magnetic helicity between 
the different scales, they will act implicitly through £ . 

From the definition of £ given above, the evolu- 
tion equation for the t.e.m.f. is d£/dt = (du/dt) x b + 
u x (dh/dt). Proceeding in a similar form as in Refs. 



l|Blackman fc FieldlliooalKandus et alj|2006l ). it now reads 



d£_ 

dt 



iv.fu + fu. 



b V x U 



(12) 



+ - [(V x b) - b-u - (V x u)J B 



--u 2 (V x B) + r/u x V 2 b 



uV 2 \i x b + T. 



where T are the triple correlations for which a closure must 
be applied. Observe that the presence of U adds two new 
terms to the equation for £ but does not explicitly modify 
the ones found in the absence of those flows. The influence 
of U on the terms proportional to B will be through the de- 
pendence of those terms w ith the magnetic helicities (cf. Ref. 
iBlackman fc Field! (2002)). To gain conceptual clearness we 
shall make further physical hypotheses on our systems, that 
will also help to simplify the mathematics. One of them is 
to consider that large scale flows U are fully helical. This is 
consistent with the concept of mean field dynamo and with 
the chosen boundary conditions. Therefore to track the evo- 
lution of the large scale velocity flow, we shall study its asso- 
ciated kinetic helicity, defined as H u — ■ , where 
W = V x U is the vorticity. The derivation of the equation 
for H u is explained in the Appendix, and the result is 



dt 



2(V x b) x b ■ W + 2^ V U ■ W. 



(13) 



where the semi-equality is due to the fact that we are ap- 
proximating volume average by a local spatial average. It is 
well known tha t kinetic helicity is not conserved in MHD 
(|Biskamdll997h , so eq. (|13p is not an essentially new result. 
However it serves to our purposes in showing that large scale 
helical flows can be induced by turbulent b-fields, provided 
they are not force-free. At this point we make another sup- 
position: we take V ■ £ — which, besides being consistent 
with the chosen boundary conditions, means that the induc- 
tion of large scale magnetic fields maximal for U = (cf. eq. 
[4]). Observe that by imposing Coulomb gauge on eq. ([5]l| 
we obtain a further constraint on the mean fields, namely 
V ■ £ = — V ■ (U x B) , and the fact that we consider it equal 
to zero allows us to replace B ■ (V x U) = U ■ (V x B) . 



3 IMPLEMENTING THE TWO SCALE 
APPROXIMATION 

As was advanced in the Introduction, we shall work within 
the two scale approximation, whereby mean fields are sup- 
posed to peak at a scale fe^ 1 , and stochastic ones at k^ 1 <§C 
k^ 1 . We begin by noting that eq. (|12|l together with the 
definition of the t.e.m.f. and eq. (|13[1 . are very compli- 
cated, involving new functions of the mean and stochas- 
tic fields for which further equations must be deduced. 
From the constraint on U and B derived from V • £ = 0, 
and the fact that we are considering B as force-free, we 
can write U • (V x B) ~ fciU • B, where the last semi- 
equality stems from the fact that B is considered to be 
force- free. From eqs. (|10p and (|lip . we see that the dot prod- 
uct of £ with B is responsible for magnetic helicity trans- 



port. Let us write £ B 



£ ■ 



Its evolution equation is 



2 We cannot use equation (8) for a to calculate (V ■ £)o as it is 
trivially zero. 

o — ~B 

J Observe that in contrast to previous works, here £ is the dot 

product of the t.e.m.f. with B, and not the component of that 
force along the mean magnetic field. We chose to work with this 
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d£ B /dt 



(d£/dt) ■ B + £ ■ dH/d t . Proceeding similar ly as 
in Refs. (|Blackman fc Field! ,20021 : iKandus et al] 120061 ), we 
obtain the following full form for the evolution equation for 

£ B , in the two scale approximation: 



9r 
at 



k L H c s H c L +\k L [klHl'-H-] \Hf 



(14) 



where we replaced U ■ B 
scale cross-helicity "' 



Hs = 



(U-B) = Hf, the large 

\ / vol 

(u ■ b) vol ~ u ■ b, the small 

u 2 /2 and 



E" 



scale cross-helicity; (V x b) • b 
H u — u • (V x u) where these two last quantities are consid- 
ered prescribed. Since we are considering B to be force-free, 
we replaced |b| 2 ~ k L \h%\ and (V x 5) • B ~ k\H% . 

The last term in eq. (|14[1 . Ci£ S j includes the effect of vis- 
cosity, resistivity, the term £ -dH/dt and, more importantly, 
the three point correlations denoted by T in eq. (|12[) . We 
see that besides the equations derived until now we need 
the ones for H^s- The derivation of these equations is 
sketched in the Appendix and, although being a straight- 
forward procedure, it is a rather tedious one and we end 
up with a system of seven equations, besides the four ones 
already shown above: the two equations for H£ s , the ones 

for F B = (V x b) x b ■ B, i.e., the dot product of the small 

— — ~w — 

scale Lorentz force with B, for F 



F U 



(V x b) x b • W, for 

(V x b) x b ■ U, for £ W = £ ■ W, i.e., the scalar 
product of the t.e.m.f. with the large scale vorticity, and for 
E b — b 2 /2, i.e., the small scale magnetic energy. Eq. (|14|) 
is already expressed in the two scale form. The other ten 
equations read 



dt 



dt 



2£ B 



-2£ 



r, 7 2 rrM 



2 V k 2 3 H^, 



dH u —W 



dt 



2F W -2vk 2 T H L 



(15) 



(16) 



(17) 



dt 



dt 



dt 



k L Hg\H u °\ + h L [k%H?-H° 



3 ' 1 3 

-2k L E u ] HE - ( 2 £ 



- F B + £ W -2{v + ri) kiHg, 



-F B - £ W -2{u + n )H<E, 




CfF 



(18) 
(19) 

(20) 
(21) 
(22) 



®F ^ 1 7 3 tjC I rjM I 

and 



-2F\ 



k s H s £ + H s F - C, b E 2 , 



(23) 



(24) 



dt 

where E% = (E b ) 



4 MAKING THE EQUATIONS 
NON-DIMENSIONAL 

In order to work with non-dimensional quantities, we de- 
fine the following dimensionless variables: r = ksut, J rW — 

F W /k%u\ T u = F^/kW, T B = F B /k s u\ E b = E b /u\ 
Q B _= £ B /u\ Q w = £ W /k S u\ Hf L = ksH^/u 2 , H~ = 
H u /k s u 2 , Hl, s = Hg s /u 2 , H u = H u /k s u 2 , E u = E u u 2 , 
£i = (i/ksu (i = 1, 2, F, b), R e = vk s /u, R m = rjk s /u, with 
R e the Reynolds number, i? m the magnetic Reynolds num- 
ber, and r = kt/ks- The non-dimensional equations then 
read 



9Q E 



dr 



duj 1 

dr 

dnf 

dr 

dt 

dU<£ 
dr 

dH c s 
dr 

dr w 

dr 

dT" r 
dr ' 

dT B 



\rH c s H c L +\r [H¥~H U ] \H 



3 13 u 3 
-r a H L -£iS , 



2 t,c 

6Q W 



H 



+ ±r [H% - H u 



2rE" 



= 2Cf-^-H¥, 

JMn 



_ on B * ~jM 

- -2Q. - — H S , 



2T V 



Re 



T-B . f^W ( 2 2 \ 2 njC 

T +Q ~ \ ~R~ + R~' r Hl > 

\ rim tie 



—T B — Q w — (— — h — 



its , 



-r U s ri L - -r yJ2=? 2 H - Cf^ 



H 



1 3t_,C \njM I 

^— - H s \H L 
dr 3 11 



and 
dr 



r 



^Hl - £ F T B 



■-2 - rt s y + rL s J- - ?6-2- 



(25) 

(26) 
(27) 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 

(35) 



5 NUMERICAL RESULTS AND DISCUSSION 



function because it facilitates the numerical integration of the 
equations, while preserving the generic behaviour of the projec- 
tions. 



We numerically integrated equations (|25[) - (|35[ ) using the fol- 
lowing parameters and initial conditions: r = 0.2, H u = 
-1, H" = 1, Q s (0) = Q w (0) = 0, Hf (0) = 0.001, 



Hf (0) 



-0.001, H (0) = 0.0001, Hi (0) 



0.0001, 
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Hg(0) = -0.0001, T B (0) = T w (0) = T u (0) = 0, 
3§ (0) = 0, Rm = -Re = 200 and 2000 (i.e., magnetic Prandtl 
number P m = 1), = 1/2 (strong non-linearities) and 
fi = 1/Rm (weak non-linearities jjj. A comment about the 
chosen values for £, is in order: In principle this parameter 
can depend on R m ; however, results of numerical simulations 
show that it is of order unity for R m < 100. In this sense, the 
value £ = 1/2 would be in accord with those results. As we 
are working here with larger values of Rm for which, to our 
knowledge, there lacks numerical estimations of £j, we chose 
two values that might represent the two extreme behaviours 
of this parameter. Nevertheless, we must stress that the va- 
lidity of this choice should be checked by direct numerical 
simulations. In Fig. [1] we plotted H u as a function of r for 
fi = 1/2. The long dashed line corresponds to R m = 200 
while the short dashed one to R m — 2000. We see that the 
generation of U is rather weak, but the effect seems to be 
stronger for R m — 2000 as time passes. In Fig. [2] we plotted 
7i u as a function of r for £j = 1/Rm. The full line corre- 
sponds to R m = 200 while the doted one to Rm = 2000. In 
this case there is a strong production of large scale kinetic 
helicity, it being stronger for R m — 2000 at the beginning of 
the integration, while for later times there seems to be no 
difference between the outcomes for the two R m considered. 
In Fig. [3] we plotted the the logarithm of the small scale 
magnetic energy, In as a function of r, for fj = 1/2. 
Each curve consists of two curves: one with the effect of U 
and the other without this field. This superposition of curves 
means that for the chosen value of the effect of U on the 
evolution of small scale magnetic energy is negligible. The 
upper curve corresponds to the largest value of R m , and we 
see that in this case a saturation value for H2 larger than 
for R m = 200 is attained. In Fig.|4]we plotted the logarithm 
of the small scale magnetic energy In (St), for fi = 1/R m . 
Long dash curves correspond to R m = 2000: upper curve 
contains the effect of U, lower oscillating curve is without 
the action of those fields. Short dashed curves correspond 
to R m = 200, with the same features for the presence and 
absence of U. We see that the action of U strongly enhances 
the generation of small scale magnetic energy, and again this 
effect is stronger for larger R m - In Fig. [5] we plotted Q B as 
a function of r for £j = 1/2. Here again each curve con- 
sists of two curves, one with the effect of U and the other 
without, showing again that for strong non-linearities the 
effect of those flows is negligible, consistently with previ- 
ous figures. Upper curve corresponds to R m = 2000 while 
lower curve to R m = 200. We see that again for larger R m 
there is an enhancement of Q B . In Fig. [6] we plotted Q B as 
a function of r for = 1/R m - We see here again that the 
action of large scale flows enhances the mean electromotive 
force Q B and this enhancement is stronger for larger Rm- 
Oscillating, dotted-line curve corresponds to R m = 2000, 
and full-line curve to R m = 200. The curves corresponding 
to the absence of the effect of those flows are almost indis- 
tinguishable from the r axis. In Fig. [7] we plotted TC^ as a 
function of r for £ = 1/2. Here again each curve consists of 
two curves, one with the effect of U and the other without, 
showing again that for strong non-linearities the effect of 

4 As before, subindex "i " denotes 1, 2, F and b. In each integra- 
tion we assume that all are the same for all "i" . 
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Figure 1. Large scale kinetic helicity H u as a function of r for 
= 1/2. The long dashed line corresponds to R m = 200, and 
the dotted line to R m = 2000. The generation of H u is stronger 
for the largest value of R m ■ 

those flows is negligible. Fast growing curve corresponds to 
Rm = 200 while lower one to R m = 2000. The coincidence 
of the two curves for short times corresponds to the kine- 
matic regime, where back-reaction of the induced magnetic 
fields b did not take place yet. In Fig. [8] we plotted Tijf as 
a function of r for £ = 1/Rm. We see here again that the 
action of large scale flows enhances the mean electromotive 
force Ttj~ and this enhancement is stronger for larger R m - 
Dashed curves correspond to R m = 2000: the ones with the 
largest amplitude correspond to the action of U, while the 
lower amplitude to the absence of this effect. Full line cor- 
respond to Rm = 200, and the features with respect to the 
presence and absence of large scale flows are the same as for 
Rm = 200. The coincidence of all four curves at the begin- 
ning of the evolution corresponds to the kinematic regime. 
In Fig. [9] we plotted Hi, as a function of r for = 1/2. 
Dashed line curve corresponds to R m = 2000 while full line 
to Rm = 200. Consistently with Fig. fTJ, we see that Hi, is 
larger for larger R m - In Fig.[l0]we plotted H*£ as a function 
of r for £j = 1/Rm- Dotted line corresponds to R m = 2000 
while full line to R m = 200. Consistently with Fig. ([2]), we 
see that Hi is larger for R m = 2000 than for R m = 200, 
with the difference in amplitudes between both quantities 
getting smaller with time. 



6 CONCLUSIONS 

In this paper we studied semi-analytically and qualitatively 
the generation of large scale flows by the action of a turbu- 
lent mean field dynamo, and the back-reaction of those flows 
on the turbulent electromotive force for two values of mag- 
netic Reynolds number, R m = 200 and 2000, and magnetic 
Prandtl number P m = 1. We considered a system in which 
small scale turbulent flows are fully helical and prescribed 
by a given external mechanism, i.e., a kinematically driven 
dynamo, and that this system possesses boundary condi- 
tions such that all total divergencies vanish. The turbulence 
was considered to be homogeneous and isotropic, which al- 
though being of limited applicability to obtain quantitative 
results for real systems, it serves to study many conceptual 
aspects of large scale magnetic field generation, besides enor- 
mously simplifying the mathematics. We followed the evo- 
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Figure 2. Large scale kinetic hclicity H" as a function of r for 
£i = l/Rm- The continuum line corresponds to R m = 200, and 
the dotted one to R m = 200. At the beginning, the induction of 
H u is stronger for the largest value of R m . 
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Figure 3. Logarith of the small scale magnetic energy Hjj as a 
function of r, for fj = 1/2. Each curve is in fact two curves, with 
and without the effect of U, which means that in this case the 
effect of those flows is negligible. The upper curve corresponds 
to R m = 2000 and the lower one to R m = 200. The small scale 
magnetic energy density is very small although larger for R m = 
2000 



lution of large scale flows through their associated kinetic 
helicity H u , as one of the suppositions we made was that 
those flows were fully helical too. One crucial assumption we 
made was that V • £ = 0, which due to the Coulomb gauge, 
imposed a further constraint on U and B that allowed us 
to expressed several terms as large scale cross-helicity. An- 
other one was to assume that B is force-free and the main 
reason to adopt it is that it helped to simplify the math- 
ematics. Although this condition can be fulfilled in certain 
astrophysical environments, this is not a general situation, 
therefore it should be dropped off in future works that aim 
at generalizing this work. 

We found that large scale flows act on the t.e.m.f. £ 
through large and small scale cross-helicities, and that, for 
the minimal r closure considered here, the effect of those 
fields is stronger for large relaxation times ( £j = 1/R m ). 
For short relaxation time (£j = 1/2), the effect of those 
fields seems to be negligible. The choice of the values for £» 
was arbitrary, in the sense that, to our present knowledge, it 



Figure 4. Logarith of the small scale magnetic energy as a 
function of r for = 1/R m - The upper, long dash curve corre- 
sponds to the presence of U, while the lower, oscillating one, to 
the absence of those flows, both of them for R m = 2000. Short 
dashed curves represents the same quantities but for R m = 200. 
In this case large scale flows strongly enhance the production of 
small scale magnetic energy, and this effect is again stronger for 
larger values of R m . 
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Figure 5. Mean electromotive force Q 8 as function of r for £j = 
1/2. Each curve is two curves, one with the effect of U and the 
other without those fields, showing that for this value of £ the 
effect of those fields is negligible. Upper curve corresponds to 
Rm = 2000 and lower curve to R m = 200, which shows that for 
the largest R m , Q is slightly stronger . 



is not known how that parameter depends on the magnetic 
Reynolds number for R m > 100. For R m < 100 it seems to 
be confirmed that £j is of order unity. Here we chose to work 
with two values that may be considered as representative of 
two extreme possibilities: & = 1/2 would be consistent with 
the predictions of the numerical simulations (although they 
were made for a different R m interval), while £ = 1/Rm 
would represent a resistive case. In any case, more reliable 
values should be given by numerical simulations performed 
for Rm > 100. 

Due to the simple system considered and the approxi- 
mations we made, we do not intended to find quantitative re- 
sults, as for example estimate the time interval during which 
U x B < £ is valid, nor do we extract more conceptual and 
qualitative conclusions. We end this work stressing the im- 
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Figure 6. Mean electromotive force Q 8 as function of r for 
= 1/Rm- Dotted line corresponds to R m = 2000 and con- 
tinuum line to Rm = 200. The curves corresponding to the ab- 
sence of large scale flows are of negligible amplitude and almost 
indistinguishable from the t axis, showing that in this case the 
enhancement of the t.e.m.f. by the shear fields is very strong. 



Figure 8. Large scale magnetic helicity as a function of r, 
for £j = 1/Rm- Dotted curves correspond to R m = 2000 while 
continuous curves to R m = 200. In each case, the strongly oscil- 
lating curves correspond to the action of large scale flows, while 
the slowly oscillations to their absence. Consistently with what 
was shown in Fig. \E\ the action of U enhances the cascade of 
magnetic helicity. 



1.2 

: 

0.8 
0.6 
. 4 
0.2 






II' 

. 025 
. 02 

. 015 
. 01 



100 125 150 175 



200 



50 75 100 125 150 175 200 



Figure 7. Large scale magnetic helicity as a function of r, 
for £i = 1/2. Again in this figure, each curve is two curves, one 
with the effect of U and the other without, showing that the effect 
of those flows on the evolution of magnetic helicity is negligible 
for the chosen value of The growing curve correspond to R m = 
200, while the slowly growing, lower curve to R m = 2000. 



portance of studying this problem via numerical simulations, 
that will show us the next paths to follow in a further an- 
alytical study, besides confirming or contesting the results 
presented here. The semi-analytical study of the anisotropic 
case is also of the most importance, as well as the consider- 
ation of other boundary conditions. 
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Figure 9. Large scale cross helicity as function of r for = 
1/2. Consistently with Fig. [T] the generation of large scale cross- 
helicity is stronger for R m = 2000 (dashed line) than for Rm = 
200 (full line). 
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A2 Evolution Equation for the Large Scale 
Kinetic Helicity 

It is obtained as dH~/dt = (u ■ dW/dt) + 
■ dU /di) . Replacing the corresponding equations we 
obtain 

8H 11 



dt 



2(V x b) x b • W - 2vk\H u , 



(A3) 



where the semi-equality stems from the fact that we approx- 



imated (■ ■ ■) j ~ ■ ■ ■. We define F = (V x b) x b-W, and 
thus eq. (|A3f) reads 



Figure 10. Large scale cross helicity as function of r, for 
£i = l/Rm- Consistently with Fig, l2lthe generation of large scale 
cross-helicity is stronger for R m = 2000 (dotted line) than for 
Rm = 200 (continuous line). 
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dt 



2F 



2vk 2 T H u . 



(A4) 



A3 Evolution Equation for F = (V x b) x b and its 
Projections 

It is found by taking curl of eq. and using it to expand 
[d (V x b) /dt] x b + (V x b) x db/dt. After a somewhat 
lengthy, but straightforward calculation, where it was as- 
sumed that for the two scale approximation V ■ F we 
obtain 



dF 

dt 



|(V-£) (VxB 
2 



ks ' 



V.S)B 



£°V^U-(>F. 



(A5) 



As in the body of the paper, we assume that V ■ £ = 0, so 
eq. (| A5|) reduces to 



APPENDIX A: DEDUCTION OF THE 
COMPLEMENTARY EVOLUTION EQUATIONS 

Here we sketch the derivation of the evolution equation for 
the large scale kinetic helicity as well as the set of extra 
equations needed to study the problem considered in this 
article. 



Al Evolution Equation for the Large Scale 
Vorticy 

We start from Navier-Stokes equation written in the form 



dt 



= -(VxU)xU-V 
+ (V x B) x B + ;A7 2 U. 



2 



(Al) 



The equation for W = V x U is obtained by simply taking 
the curl of eq. (|A1|) . after replacing the decomposition in 
mean and stochastic fields, and using the hypothesis that u 
is fully helical. We have 



aw 

dt 



-V x (W x U) + V x (V x b) x b 
uV 2 W. 



(A2) 



dF 
~dt 



3 b 



-E b X7 2 \J- 

3 



CfF. 



(A6) 



To find the evolution equations for the scalar product 
of F with W and B, we use the above defined expression 

F , and an analogous expression for B. Again the evolu- 
tion equation is found by taking the time derivative of the 
complete expression. In the two scale approximation we have 



dF W 
dt 



-^HsV 2 B ■ W + ^E b V 2 U ■ W - ( F F W 



tfHgB ■ W - ^E b H u - CfF W . 



(A7) 



V • (U x B) 



0, we can 

c 



Due to the fact that V • £ = 

write B • W ~ U ■ (V x B) ~ k L U ■ B ~ kt,H%. Using 
the fact that for a fully helical U field we can write | W | ~ 

1/2 



1/2 



A: 



dF 



we obtain 



dt 



lh± E b H u 
3 



(fF 



(A8) 



For the projection of F along B and along U we proceed 
— w 

analogously as for F . Using the fact that for a large scale 
force-free field we can write IBqI = kV 2 l-H^ 7 ! 1 ^ 2 , we obtain 



V ■ ((V x V) x 6) 



-b- V 2 6- IV x b\- 



:k 2 s \b\ 2 



k 2 \b\ 2 
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9F k L a I rjM I 
— _ -H S \H L | 

and 



dF U kl „ c „c 2k L b 

— _ IT q IT j — Hi 

dt 3 s 3 



where we used U ~ \H U \ /k L 



(A9) 



(A10) 



between spatial averages of functions of stochastic and mean 
fields. 

dE b 



. )± - -(V x b) x b ■ U+ (V x b) x u B - ( b E b . (A17) 

To deal with the second term we write V x b = 
[(V x b) • b] b/ |b| 2 - [(V x b) x b] x b/ |b| 2 and thus 



(V x b) x u 



[(V x b) ■ b] b x u 



A4 Evolution Equation for the Cross-Helicity 

Cross-Helicity is defined as H c = (U-B) tl0 ;. After obtaining 
from eq. (|A1|| the evolution equations for U and u and using 
eq. @ and Q, we obtain the following equation for the large 
scale cross helicity, Hjj and the small scale cross helicity : 



and 
dl Hs 



(Vxb)xb-B + £:-W 

2 (v + V ) k 2 L HZ 



-(Vxb)xb-B-£-W 

2(^ + 7?) klHt 



(All) 



(A12) 



+ 



{[(V x b) x b] x b} x u 



(Vxb)-b£ {[(V x b) x b] ■ u}b 



(V x b) x b u • b 



We obtain for the second term in eq. (|A17|I : 



(V x b) x u ■ B 



+ 



(V x b) -b£-B 



u b (V x b) x b • B 



(A18) 



Replacing F B = (V x b) x b ■ B, defining £ = £ ■ W, and 
using the fact that for fully helical U we can write |w| ~ 

1/2 

H u , in the two scale approximation we have 



k 

dH L 
dt 
and 



dt 



£ +F -2{v + rj)kiHZ 



-£ W - F B - 2 (y + n ) k 2 s H%. 



(A13) 



(A14) 



k 2 s H^£ a Hg~F , Ain . 

where the second term in the second row of expr. (j A18f) 
was considered to give a null contribution when averaged. 
Defining E% = (E b ) we can write the evolution equation 
for the small scale magnetic energy as 



dt 



-2F\ 



/ E b 



-klHf£ B + H c s F B 



(A20) 



A5 Evolution Equation for £ = £ ■ W 

Using equation (JX3J , and V • £ — 0, we obtain 



— w 

i \h c s |w| 2 + ^(fcitff -h«)h£ 



dt 



i« 2 (V x B) ■ W - ( 2 £ l 



(A15) 



where in the last term we considered the term £ ■ dW/dt, 
and the three point correlations. Performing ( V x B) ■ W ~ 
k L B ■ W ~ k L (V x B) ■ U ~ k\H% , where the semi- 
equality before the last stems from the fact that V • £ = 
- V ■ (U x B) ~ 0, we obtain 



dT 



dt 



2 r , 



t't, I 



+ \k L 



— w 

- (2£ . 



[k 2 s Hf -H u - 2k L E u ] Hi 
(A16) 



A6 Evolution Equation for E b 

It is obtained by scalar multiplying eq. Q by b and then 
taking volume average. In order to simplify the mathemat- 
ics, we approximate the volume averages by a dot product 



